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Abstract 

This paper introduces a new difference scheme to the difference 
equations for N-body type problems. To find the non-collision periodic 
solutions and generalized periodic solutions in multi-radial symmetric 
constraint for the N-body type difference equations, the variational 
approach and the method of minimizing the Lagrangian action are 
adopted and the strong force condition is considered correspondingly, 
which is an efficient method in studying those with singular potentials. 
And the difference equation can also be taken into consideration of 
other periodic solutions with symmetric or choreographic constraint 
in further studies. 

1 Introduction 

The N-body problem is a classical and important problem in celestial me- 
chanics and mathematics. It consists of determining the orbits of N bodies 
interacting in accordance with the gravitational law of Newton. That is two 
bodies attracting each other, the force of attraction being directed along the 
line joining them, proportional to the product of the masses, and inversely 
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proportional to the square of the distance between them. For given N bodies, 
the preceding N-body problem can be described by the following nonlinear 
system of second order differential equations: 



A' 



mm = 

l=X,l^i 



I = 



1,... ,N, 



(1) 



where rrii is the mass of zth body, qi is the position of zth body in R d , N is 
the number of bodies, and G is the universal gravitational constant, which is 
always taken to be 1 for convenience in mathematics. Hence we take G — 1 
in this paper. 

Poincare showed that it is impossible for N > 3 to find an explicit ex- 
pression for the general solution. So even for the three-body problem, it 
is impossible to describe all the solutions. Till now, for planar three-body 
problem, there are only three orbits are found clearly for three equal masses, 
including the existence and the shape of the orbits. They are Euler's solu- 
tion (see FigH]), Lagrange's solution (see Figj2]) and Figure-Eight solution 
(see Figj3]). The Euler's solution was discovered by Euler in 1765. It de- 
scribed a collinear configuration, in which one body locates on the center of 
masses of the three bodies, and the other two rotate around the center of 
the masses on the two endpoints of the diameter of the circle orbit. The La- 
grange solution was discovered by Lagrange in 1772. The three bodies form 
an equilateral triangle on a circle orbit and the center of masses is the origin 
all the time. The Figure-Eight solution is a new, surprisingly and remarkable 
solution after the two solutions, which was found numerically by Moore first 
in 1993 (see [5]). In 2000, Chenciner and Montgomery rediscovered it and 
proved its existence at the same time. It describes three equal masses chas- 
ing one by one on an Eight-shaped curve in the plane. We refer the reader 
to the following articles and references therein for more detailed description: 
Chenciner et al [2], Montgomery [6] and Simo [8] etc. 

There is no other specific result for N > 3. Hence the numerical solution 
is more and more important in finding and describing new solutions vividly. 
In recent years, some new periodic solutions of the N-body problem have 
been found by minimizing the action of cycles with some special symmetric 
or choreographic constraint (e.g. [3] , [TTJ etc), in which Fourier approach is 
used by most authors for the numerical solutions of N-body problem. In this 
paper, we propose the N-body type difference equations, and from which 
we can discover some new periodic orbits by minimizing the corresponding 
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action functional. 





Figure 1: Euler's Solution 



Figure 2: Lagrange's Solution 



Figure 3: Figure-Eight Solu- 
tion 



The N-body type problems have been studied by Bahri, Rabinowitz [I], 
Zelati [TD] among others. In [10j, Zelati proved the existence of periodic 
solutions under the radial symmetry assumption. In this paper we prove the 
existence of periodic solutions for the difference equations with multi-radial 
symmetric constraint, which includes the radial symmetry. In the section on 
numerical examples, we propose a new variational difference method for the 
periodic solutions numerically. 

The periodic solutions for N-body type problems (including N-body prob- 
lems) can be described by 



where N, rrii and q^ are defined as above, q = (qi, - ■ ■ , ?at) in (R d ) N , T is the 
period, and generally the singular potential function U is 



miqi = VU(q), i = 1, • - • , N, 
q(0)=q(T), q (0) = q(T), 



(N) 
(P) 




(Al) 



oo, as |£| — > 



(A2) 



U(q u --- ,q N )<0 
where Uu is usually described by 





vtiimi 



(2) 



Qi-Qi 



a 
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2 Difference Equations 



In order to formulate the difference equations for N-body type problems (1N|) - 
([P]), we discretize time within one period [0,T] with k + 1 points to, • ' ' 
and denote the set of these discretized points as S, i.e. S = {t , • ■ ■ , tk}, and 
introduce a difference scheme by 



q(t) = q(t) :-- 



-q(s-l) + - — T^qis) 



h 



h 



k 



(DS) 



where q(s) is an abbreviation of q(t s ), t s is the sth discretized point of [0, T], 
Q — (^i, ■ • ■ ,?at), and h is the discretized step. Moreover, the traditional 
forth difference scheme is used to approximate the derivative of q, which is 
used in [9]. Then the first order difference of q is 

Aq(t s ) = q(s + l)-q(s), 

and the second order difference of q is 

A 2 q(t s ) = q{s + 1) - 2q(s) + q{s - 1). 

Hence the derivatives of q(t s ) can be approximated by 

q(t s ) = -Aq(t s ), q{t.) = ^A 2 q(t s ), 



By the difference approach above, (TNJ) can be transformed into 
—^rriiiq^s + 1) - 2&(s) + q^s - 1)) = V q . {s) U, 
where i — 1, • • • , N and 

U = Y J U{q{s))= l -Y. E tf«(«oo -*(*))■ 

For convenience, we set the following notations: 

/ 2 -1 ■■■ -1 \ 

-l 2 -l '•. ; 

M k = o ' • • ' • • ' • • 



V 



-1 2 -1 
0-12/ 



(3) 



fqj\ 



( ^ \ 



i(2) n 



and 



M 



where 



Mi 



( Mi 

V 

/ rriiMk 



(d) 1 



I a ip)T \ 

?i(2) 



\ o (p)T / 



M N J 



h 



N, p 



rriiMk 



It is easy to verify that M is a positive definitive matrix. Then the difference 
equation of (jNj) becomes 



-M • g - V Us) Uh = 0, 



(Nl) 



and it is obviously that the difference boundary conditions can be described 
by 

q(j) = q(k + j), ieZ+, (Pi) 

where Z + is the set of positive integer. Hence the Lagrangian functional of 
l|NIjl - (|FIj» is the following 



J(q) 



2h 



< Mq,q> -Uh. 



(4) 



Theorem 1. There exists a small number S > for alii ^ j, i,j = 1, • • • , N, 
s G S, such that \qi(s) — qj(s)\ > 5. Then the solution q of $Nl\) - [Pl\) 
approximates the solution q of $M)- $B) when k — > oo. 



Proof. Since \qi — cjj\ > 5, we can deduce the following equation from (|2J) and 
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1 k 

?E E 



s=l \<i^l<N 
k 



l9!W-?iWI (o+1) 

< iE E 

= c, 



(5) 



where C is a constant. Similar to (IDSj) . we introduce the difference scheme 
to i.e. we have 



-q(t.-i) + 



h h 
where t G [t s _i, t s ] and t s _i,t s G 5. Notice that 



q(t s 



?(*) = ^(?(*«) - ?(*-0) = ^Ag(t,_i). 



Then when we let 

x = (g,g), 

the system (JHD-© 

can be rewritten as follows: 
x = F(x), x(0) = x(T), 
where x satisfies (JDS) and F G C(R 2dN , R 2dN ), that is 

/ ; \ 



-£Vl7( ff ) 



V 



/ 



From (jSJ), we have |F| < Mo for some constant Mq, and for Ve > 0, there 
exist 5' > 0, such that 

|F(x!)-F(x 2 )| <e, 



6 



when \xi — x 2 \ < 5' '. From (IDSp . for t e [t a -i,i a ], s = 1, • • • , k, clearly we 
have 



\xk{t) - x(t a -i)\ 



t s -t 



h 

t — t s _l 



x(t s -i) + - — ^-x(t s ) - x(t s -i) 



h 

t — t s -\ 



h 



(X(t s ) - X(tg-l)) 

F(x(t s ^))(t s -t s ^) 



(6) 



= |(* — t s _i)^(^(t s _!))| 

< hM . 

Therefore, there exists a K(e) > such that 

\F(x k (t)) - F(x(t s _i))\ < s, whenever k > K(e). 
The error within one period can be estimated by 

A fc (t) = x k (t) - (x(0) + jf F(x fc (r))rfr 

Hence 

k „i 

A *( T ) = / ' (F(x(ts-i))-F(x k (r)))dT. 

S=l Jta-1 

From (Q, noting that the discretized step is h — 1, we have 

k 

\A k (T)\<J2^ = Te. 

s=l 

This means A k — > 0, — > x, i.e. q —>■ q when — > oo. The theorem is 
proved. □ 

Let 

Q = {q\qe (R d ) N , s.t. q t ? q u Vi ^ Z}, A = {<? | <? G H l {R/T\ Q)}, R d = JJi^ 



i=l 



Then the multi-radial symmetry can be described as follows: 



Ao = {<? | <? e A, <f l) (s + = z = 1, • • • , b d , 6< = d}, (Ms) 



where k % = A/A4, A4, A G Z + , A is a common multiple of A4, and q( bi '(s + 
k l ) = —q( bi \s) means 



q^)(t) = -q^\r), where t G [t,^, t a+%i \, r G [t s _ u 



3 The Existence of Periodic Solutions 

Let H 1 is a Hilbert space, q G -f/' 1 ([0, T], (R d ) N ), i.e. ft G i^flp, T], R d ), 
i — 1, • • • , N, and we define the norm of q in H 1 by 



II q ||= V 7 ^ 2 + S 2 , 

where 

^=^e I s \^)\dt, s 2 =j2 r \m\ 2 dt. 

s=l Jts-1 s= l Jt B -i 

The following lemma describes the relationship between the critical points 
of J and the periodic solutions of fINip - flPlI) . 

Lemma 1. (g(l),--- , q(k)) T is a critical point of J if and only if (q{0) , ■ ■ ■ ,q(k+ 
1)) T is a solution of Wl\)-lPl\). 

Proof. Since the difference equations transform the infinite dimensional dy- 
namics into finite dimensional dynamics, it is reasonable to view J as a 
continuously different iable functional on H 1 . 

Let q = (5(1), ■ • • , q(k)) T is the critical point of J if and only if the first 
order derivative of J on q is equal to zero. That is 

f(q) =Mq- V 4 U = 0. 

Notice that 

(g(0), q(l), • • • , q(k), q(k + 1)) T = (q(k), q(l), • • • , q(k), g(l)) T . 
This satisfies ( INlft - llPip obviously. Thus the proof is completed. □ 

Remark 1. Actually, is a k-polygon in R d as a simple case with the 
difference scheme liDS\) . Hence we can call {<?( s )}sio = (<?(0), ■ • • , q(k + 1)) T 
a 'cycle' in (R d ) N , which means obviously that {(^(s)}^ = (^(0), • ■ ■ , qi(k + 
1)) T , i — 1, • • • ,N, are 'cycles' in R d . 
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We denote 

k 

A = (J{t G [t s _i,t s ] I qi(t) = q~i(t), for some i ^ I}. 

s=l 

Definition 1. We call q is a non-collision solution of liNl\) and liPl\) . if 
A = 0. 

Definition 2. We call q is a generalized solution of liNl\) and liPl\) . if 

a) meas(A) = 0; 

b) q satisfy / fMl) : 
c) 

~ U(q(t)) = C, (7) 
i=i 

where t e [0,T]\A and C is a constant. 

Remark 2. It is easy to verify that the integral of ([?|) with respect to t 
approximates J(q) since 

/ U u (qi(t) - Ut))dt = U u (qi(s) - ft(s)). (8) 

For convenience, we suppose U in (1N1I) has only one singular point, and 
without loss of generality, let it be G R d . Hence we have 

Definition 3. Let some V a G ^(^{O}, R), V a +oo (£ -> 0) and 

-*/«(£)> I W«(fl| a , VeGi? d \{0}, (SF) 
where \£\ is small. Then we say Uu satisfies Strong Force (SF) condition. 

The condition (ISFI) is very efficient to overcome the obstacle coming from 
singular potential. In [3], (ISFI) is used to bound q away from S, where S is 
a closed nonempty set, and U — > — oo as q — > S. Here, we introduce Lemma 
2 based on the discussions above for N-body type difference equations. 

Let T be the family of homotopic 'cycles' in H X (%T\, (^{0})^), then 
we have 
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Lemma 2. Assume condition liSF]) is satisfied, and there exists a constants 
c such that 

(Mq,q)<c, \U\<c, 

for all q G T, then q is bounded away from {0}. That is, there exists 5 > 
such that no cycles in T intersect the 5 -neighborhood of {0}. 

Proof. From (1SF|) . it is easy to verify that Uu are all negative, then from Q 
we have the following estimate for U: 

k 

\ & \ = £E E \Uu(Qi(s) - qi(s))\ < c. 



8=1 l<i^«<AT 



Let 



/ mi 



5 



\ 



mi 



\ m N J 

rh = min{mi, • • • , m N }, y = Bq, 

and let 

G = {p= (y, z) G (R d ) N x R | MG (^{0})^, 2 = V(q)} 
be the graph of V in (ISFj) . Suppose the binary relation is 

P={y,V{q)), 



10 



whose arc length is a = a(p). Then we have 

T k t 

a = J Q \d\dt = £ \d\dt 

s=l 

= t ft. VWTW^Wdt 

s=l 

< ZflAlyl' + ieivvffdt 

s=l ' ' 

< tit^mHi + ^wvi^dt 

s=l 

< E (/£ 1 H 2 ^) 5 (j£ , (1 + M VV \ 2 ) dt ) * from H51der inequality 

s=l ^ ^ ' 

1 

< (e h 2 ^) 2 (e jj_ x (i + m^v?) ^ 

= (e it, i^r*) 2 ( T + ^ £ (l w|2) rft ) 2 

< (E l^| 2 eft) 2 (t + & E J?_ a 1^) 2 from (EE) condition 
= (t\B(q(s)-q(s-l)ry (t + £ £j£jl7 ■($(*)) |dfY 

/ k N \ 2 

= EE^I(M-^-i))l 2 

' k 

*=1 l<i^l<N 



= ^ Mg~l)^T + ^\U\ 

< y/c^T + 4a c = C (a constant). 

Now let {q u }^ =1 be a sequence in I\ If q u are not bounded away from 
{0}, for We > 0, the e-neighborhood of {0} in R d , B e , must intersect with 
infinite many q v . By passing to a subsequence, we can suppose that 

inf|<f| — > 0, as n — > oo. 

This implies that the corresponding p v ascend infinitely far up to the skylight 
of G at {0} as n — > oo. 
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Suppose there exists 5 > 0, B$ is (^-neighborhood of {0}, such that por- 
tions of q u fall outside the ball B$ and the others (infinite many) q u in B$. 
Then the variation of V(q u ) would become infinite as n — > oo, but a(p) is 
greater than the variation of V(q u ), a contradiction to the fact a < C . 

On the other hand, if there is no such 5, q u will collapse to the point {0}, 
which implies that 



This contradicts to the hypotheses at the beginning. Hence U is bounded 
away from {0}. 

Now, suppose q u converge uniformly to q* in the norm in H . By the 
above argument, if q v are not bounded away from {0}, q* will be attached 



Remark 3. 

(a) For the case that the potential U has more than one singular points 
has the same result in Lemma 2. 

(b) Lemma 2 means that there is no q G V intersects the 5 -neighborhood 
of {0} under liSF]) condition. 

(c) The converse-negative proposition of Lemma 2 shows that if q e <9A, 
qt* — > q weakly in H 1 and strongly in C° , then J(g M ) — > +oo. 

With respect to the periodic solutions of (INI I) - (IP lj) . we have the following 
theorem: 

Theorem 2. Let Uu and U in (Ejj satisfy (fAl) ) -(H3] ) and q satisfy \Ms\) . 
Then j{Nl\) - fPl\) has infinitely many generalized solutions; furthermore, if 
ASF]) holds, then IIN1]) - [P1\) has infinitely many non-collision solutions. 

Proof. It is obvious that Ao is a subset of A, and U is an even function on 
A for Uu = Uu, so we can easily check that the critical points of (jSJ) on A 
are also the critical points on A. 

Since Uu does not satisfy fISFp . we can modify it as 



U 



+oo. 



to {0}. The proof is complete. 



□ 




I A A IT ) 

\Qi ~ Qi\ 



where 5 is small and 




when \qi - q { \ < -, 
when \qi — cji\ > 5. 
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Then we get a new Lagrangian functional with condition (ISFj) from (j3J), 

J 5 (g) = ^<Mg,g>-f/ 5 , (9) 

where 

& s = t E ^ 

3=1 l<i<l<N 

Let J* 5 have infimum 

cj = inf{ J 5 (g) | g G A }. 

Consider minimizing sequence g**(s) G Ao such that J 5 (g M (s)) — > q when 
/Lt — >■ +00. Then, when /x is large enough, we have 

■/ 5 (<f(s)) < c 5 + l. 
Since g G A , we have that 

< dxmax{|g^)|,|gf( S )|} 

= d x max{l|gf( S ) - #( a + ~|gf ( a ) - gf {s + fc)|} 

< d x max{-^ |g 4 (a - j) - &(s - j + 1)|, 

j=0 



^ fc-i 

-^^(s-^-g^s-j + l)!} 



j=0 

fc-i 



, fc-i 

j=0 
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where m = min |m<|. This illuminates that there exists a q G An and a 

i=l,...,jV 

subsequence {g^}^ such that 

g M -> g(s), s = 1,--- ,k, 
J\q) = eg. 

Then q is a critical point of (J3]). 

Now, from the remarks of Lemma [2] and Lemma [TJ there exist a sub- 
sequence {g M (s)} and g| in ^(R/T, (R ) ), such that g| is a non-collision 
periodic solution of N-body problem (1Nljl - (1Pljl . 

Along the line of Theoreml.l in [10], we can find a generalized solution 
q* of (|Nip - (]Pip by constructing compact nested interval cover of [0, T]\A to 
show II s — > [/ and g| — > g* when 5^0. We can also prove Qt/c is a 
periodic solution of ( 1N1[) -( |PT|) if Q is, where C is a positive integer, then we 
can find infinitely many T-periodic solutions of ( 1Nll) -( jPTI) . This completes 
the proof. □ 



4 Numerical Examples for the N-Body Prob- 
lem 

In this section, we present some numerical solutions of N-body problem ([1])- 
(TP]) as an example of N-body type. In this case, a = 1 in (j2J), which is usually 
called Newtonian Potential. 
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The method of steepest descent is adopted in finding the minimum of the 
functional J since this method is simple, easy to apply, and each iteration 
is fast. If the minimum points exist, the method guarantees that we can 
locate them after at least an infinite number of iterations. We implement 
it in MATLAB7.0 for getting the following orbits. And from the process 
in deducing J, we named the method variational difference method. The 
following examples show that it is at least very valid for our problem. 

When we set the multi-radial symmetry with 




which is a simple radial symmetric constraint in the plane. We get the 
following two figures with different initial values for three-body problems 
with equal masses. 

In the following tables, qi in the columns is the initial value of zth body, 
where % = 1,2,3. And the jth raw is corresponding to the initial value of 
the jth difference point, where j = 1, • • • , k. We set the values at random 
by MATLAB. 





<?2 


<?3 


Qi 


Q2 


93 


4.8375 
4.3081 


-3.5004 
-5.8243 


-1.3367 
1.5174 


4.8375 
4.3081 


-3.5004 
-5.8243 


-1.3367 
1.5174 


5.0723 
3.9921 


-3.3872 
-5.7708 


-1.6835 
1.7809 


5.0723 
3.9921 


-3.3872 
-5.7708 


-1.6835 
1.7809 














-4.2971 
-4.8963 


3.7003 
5.8687 


0.5988 
-0.9714 


-4.2971 
-4.8963 


3.7003 
5.8687 


0.5988 
-0.9714 


-4.5792 
-4.6106 


3.6051 
5.8569 


0.9749 
-1.2464 


-4.5792 
-4.6106 


3.6051 
5.8569 


0.9749 
-1.2464 



Table 1: The initial value of Figg] Table 2: The initial value of Fig[5] 
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Figure 4: m = [1, 1,1], J = 0.493 

When the distribution of the masses 
different numerical solutions are 
three-body problem. 





<?2 


<?3 


0.9977 
0.3897 


0.7023 
0.5013 


0.9878 
0.5027 


0.9878 
0.5027 


0.3214 
0.2917 


0.8118 
0.2859 








0.8906 
0.4815 


0.0109 
0.0134 


0.1405 
0.1421 


0.5134 
0.4002 


0.0484 
0.0046 


0.1205 
0.3442 



Table 3: The initial value of Fig. \E\ 



Figure 5: m = [1, 1, 1], J = 0.512 

are mi = 1, m 2 = 0.5, m 3 = 1.5, four 
with different initial values of J for 



Qi 


<?2 


<?3 


0.2916 
0.3066 


0.2886 
0.3100 


0.5938 
0.3580 


0.0974 
0.7207 


0.9711 
0.1348 


0.1567 
0.9382 








0.6674 
0.1623 


0.6226 
0.6986 


0.3063 
0.3450 


0.3109 
0.0311 


0.1326 
0.8893 


0.6609 
0.8840 



Table 4: The initial value of Fig. 
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Figure 6: m = [1, 0.5, 1.5], J = 0.805 





<?2 


<?3 


0.5499 
0.6606 


0.8989 
0.8988 


0.7706 
0.2469 


0.0424 
0.2721 


0.0684 
0.7092 


0.6264 
0.1738 








0.7811 
0.7249 


0.2253 
0.9763 


0.1856 
0.1436 


0.8518 
0.3835 


0.2435 
0.5769 


0.2441 
0.7770 



Table 5: The initial value of Fig. \E\ 




-6 -4 -2 2 4 6 



Figure 8: m = [1, 0.5, 1.5], J = 0.795 



Figure 7: m = [1,0.5, 1.5], J = 0.670 



Qi 


<?2 


?3 


0.9342 
0.4225 


0.3759 
0.6273 


0.1146 
0.4319 


0.2644 
0.8560 


0.0099 
0.6991 


0.6649 
0.6343 








0.1556 
0.6252 


0.6555 
0.7036 


0.7505 
0.6802 


0.1911 
0.7334 


0.3919 
0.4850 


0.7400 
0.0534 



Table 6: The initial value of Fig. E] 




-4-3-2-101234 



Figure 9: m = [1,0.5, 1.5], J = 0.832 
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The following orbits are also found for the planar N-body problems, where 
m shows the number and the masses of the bodies, J is the corresponding 
value of Lagrangian functional action. 




-4-3-2-1 1 2 3 4 

Figure 14: m = [1, 0.5, 1, 0.5], J = 1.086 



-6-4 -2 2 4 6 

Figure 15: m = [1,0.5, 1,0.5, 1], J = 1.827 



For spacial N-body problems, 



g (i,2,3) {t + ^\= -q(t)^\ E i? 3 , 



i.e. Ai = A 2 = A 3 = 2 in (IMsj) . it is also a simple radial symmetric constraint 
but in space. We find the orbits in the following figures. In each figure, the 
picture on the top left corner is the projection of the orbits on x-y plane, the 
picture on the top right corner is the projection of the orbits on y-z plane, 
the picture on the bottom left corner is the projection of the orbits on z-x 
plane, and the picture on the bottom right corner is the orbits in space. 











Figure 16: m = [1, 1, 1], J = 0.783 



Figure 17: m = [1, 0.5, 1.5], J = 0.832 
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Figure 22: m = [1, 1, 1, 1, 1], J = 2.596 Figure 23: m = [1, 2, 4, 8, 16], J = 43.348 



When we set Ax = 4, A 2 = 2 in ([Ms]) . R d = R 3 = R 2 x R 1 , we have the 



following multi-radial symmetric constraint 



«P(t + f) = -#'(*), ft (3| etf. 



In this case, the following orbits are found. 




21 



Figure 26: m = [1,2, 11], J = 5.716 



Figure 27: m = [1, 1, 10], J = 3.216 



5 Conclusion 

The N-body type difference equation is first introduced in this paper with 
a new difference scheme. To prove the existence of periodic solutions for 
N-body type problem, the variational approach and the method of minimiz- 
ing the Lagrangian functional action, both of which are typical in studying 
N-body dynamical systems, are employed in this paper. In finding the gen- 
eralized solutions and non-collision solutions numerically, we use the method 
of minimizing the functional action of N-body type difference equation. The 
presented numerical solutions of N-body problem illuminate that the vari- 
ational difference method, which is combined with the traditional steepest 
method, is rather effective in finding the periodic solutions for N-body type 
problem numerically, with these numerical examples, we can conclude that 
there exist some interesting orbits under multi-radial symmetric constraint; 
and surprisingly, even under radial symmetric constraint we have obtained 
some interesting orbits. Furthermore, difference equations transform infinity 
dimension dynamical systems into finite dimension dynamical systems, and it 
can be used in the studies for other solutions under corresponding symmetric 
or choreographic constraints. 

References 

[1] Bahri, A. and Rabinowitz, P. H.: Periodic solutions of Hamiltonian sys- 
tems of 3-body type. Ann. Inst. H. Poincare Anal. Non Lineaire Vol 8, 



22 



561-649 (1991) 



[2] Chenciner, A., Gerver, J., Montgomery, R. and Simo, C: Simple choreo- 
graphic motions of N bodies : a preliminary study. Geometry, Mechanics 
and Dynamics 287-308 (2002) 

[3] Ferrario, D. L. and Terracini, S.: On the existence of collisionless equiv- 
ariant minimizers for the classical n-body problem. Invent, math. Vol 
155, 305-362 (2004) 

[4] Gordon, W. B.: Conservative dynamical systems involving strong forces. 
Tran. A. M. S. Vol 204, 113-135 (1975) 

[5] Moore, C: Braids in classical gravity. Physical Review Letters Vol 70, 
3675-3679 (1993) 

[6] Montgomery, R.: A new solution to the three-body problem. Notices Am. 
Math. Soc. Vol 48, 471-81 (2001) 

[7] Simo, C: New families of solutions in N-body problems. European 
Congress of Mathematics Vol I, 101-115 (2001) 

[8] Simo, C: Periodic orbits of the planar N-body probelm with equal masses 
and all bodies on the same path. NATO Advanced Study Institute (Bris- 
tol: Institute of Physics Publishing) 265-284 (2001) 

[9] Yu, J., Guo, Z. and Zhou, X.: Periodic solutions of second order self- 
adjoint difference equations. J. London Math. Soc. Vol 71, 146-160 (2005) 

[10] Zelati, V. C: Periodic solutions for N-body type problems. Annales de 
l'lnstitut Henri Poincare-Analyse non lineaire, Vol 7, 477-492 (1990) 

[11] Zhang, Shiqing, Zhou, Qing and Liu, Yurong: New periodic solutions 
for 3-body problems. Celestial Mech. Dynam. Astronom. Vol 88, 365-378 
(2004) 



23 



